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Mollow triplet of a quantum dot near a metal nanoparticle 
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We present a theoretical study of tlie resonance fluorescence spectra of an optically driven quantum 
dot placed near a single metal nanoparticle. The metallic reservoir coupling is calculated for an 8-nm 
metal nanoparticle using a time-convolutionless master equation approach where the exact photon 
reservoir function is included using Green function theory. By exciting the system coherently near 
the nanoparticle dipole mode, we show that the driven Mollow spectrum becomes highly asymmetric 
due to internal coupling effects with higher-order plasmons. We also highlight regimes of resonance 
squeezing and broadening as well as spectral reshaping through light propagation. Our master 
equation technique can be applied to any arbitrary material system, including lossy inhomogeneous 
structures, where mode expansion techniques are known to break down. 

PACS numbers: 42.50.Pq, 78.67.Bf, 73.20.Mf 



Ph. 



> 

in 

o 

o 
en 



X 



I. INTRODUCTION 

The study of quantum light-matter interactions near 
metals can be used to explore fundamental quantum op- 
tical regimes such as modified spontaneous emissiort^ and 
the strong coupling regime^"—, with applications rang- 
ing from single-photon transistors^ to plasmon lasing and 
spacing^"—. Metal structures enable surface plasmon po- 
laritons which give rise to pronounced resonances in a 
similar way to high-Q (quality factor) cavity structures. 
However, metals are significantly more complicated to 
model because of material losses, and, e.g., standard 
mode expansion techniques that are well used in quan- 
tum optics theory typically fail. This motivates the need 
for quantum optics models that can be applied to metal- 
lic environments, and can expect to access new excitation 
regimes that are unique to plasmonic systems. 

Recently, there has been interest in the coherent excita- 
tion of a single atom or quantum dot (QD) near a metal 
surface. As is well known from atomic optics, a reso- 
nantly driven atom (or QD) can yield a "Mollow triplet" 
for the incoherent spectrum if the coherent Rabi oscil- 
lations have a frequency that is larger than the decay 
rates in the systeni^iiS. Ridolfo et ali^ modelled QD 
metal-nanoparticle (MNP) interactions through a master 
equation (ME) approach by assuming a single Lorentzian 
response for the metal and estimated the dot-metal cou- 
pling parameters from electromagnetic simulations; such 
an approach is useful but restricted since it essentially 
ignores the higher-order plasmon modes and cannot be 
used if the dot is too close to the metal surface — typically 
restricted to separation distances greater than a radius of 
the particle^ or else coupling to higher-order plasmons 
becomes important. Gonzalez-Tudela et ali^ employed a 
time-convolutionless (i.e., time local) ME approach and 
explored the coupling for a driven QD near a metal pla- 



nar surface; this latter method allows one to incorporate 
the full non-Lorcntzian lineshape of the metal reservoir 
for QDs that are close to the surface, however the ef- 
fects of internal couplingi^ and spectral reshaping due to 
QD/MNP-to-detector propagation were neglected. In- 
ternal coupling refers to coupling to the photon reservoir 
at the dressed-state resonances, since in general the scat- 
tering rates and radiative coupling depend on the driving 
fieldi^ii^. For a suitable MNP environment, the range 
of energy shifts of the dressed states can be substantial 
compared to the range of energies over which the LDOS 
varies. 

In this work, we introduce a powerful ME technique 
that allows one to model the quantum light-matter in- 
teractions for any general photonic reservoir function, 
including lossy inhomogeneous structures. The only re- 
striction we use is the second-order Born approximation, 
which is valid for the weak QD-plasmon coupling regime 
that we consider. We apply this approach to study the 
incoherent spectrum that is detected when a QD is driven 




FIG. 1: (Color online) Schematic showing a pump field excit- 
ing the QD metal-nanoparticle system (not to scale) , resulting 
in resonance fluorescence that can be detected in the far fleld. 
The radius of the MNP is a = 8 nm. The distance from the 
MNP to the detector is varied in the results presented below. 



resonantly near an 8-nni MNP. A schematic of this ex- 
citation scheme is shown in Fig. [1] In the strong-field 
excitation regime, we compute the fluorescence spectrum 
at a detector remotely located from the driven QD-MNP 
system, fully taking into account the effects of light prop- 
agation and optical quenching. We demonstrate that, as 
the field strength of the drive is increased, the ensuing 
MoUow triplets becomes highly asymmetric. The Mol- 
low triplets gives direct access to the regime of quantum 
nanoplamonics and contains signatures of the MNP's 
photon bath function. 

The paper is organized as follows. In Sec.|ll]we intro- 
duce the theory and ME technique for modelling a coher- 
ently driven QD in the vicinity of a MNP. We also present 
an expression for the incoherent spectrum in terms of the 
medium Green functions which arc computed exactly. In 
Sec. mil we present calculations of the Green functions 
and incoherent spectra for various pump intensities and 
pump laser detunings. The ensuing Mollow spectrum 
is seen to be highly asymmetric and we show the im- 
portance of including internal coupling effects. We also 
study the effect of the QD-MNP separation and observe 
squeezing and anti-squeczing of the spectral resonances. 
Theoretical expressions for the linewidths in terms of 
the LDOS help to explain the physics. Conclusions are 
offered in Sec. IIVI We also include an Appendix that 
presents the optical Bloch equations for this system and 
discuss useful analytical limits of the incoherent spec- 
trum. 



II. THEORY 

A. Photon Green Function 

We first introduce the classical photonic Green func- 
tion of the MNP, where the nanoparticle is assumed to 
be in air. Defining the MNP complex permittivity as 
eMNp(r,w) = e_R(r,w)+ie/(r, w), then the photon Green 
function satisfies the follow equation {uj is implicit): 

V X V X G(r,r') -fc2£MNp(r)G(r,r') = kl5{r-v')l, 

(1) 

where fco = w/c, /i = 1, and I is the unit dyadic. We 
use a Drude model for a silver MNP, where £{uj) = 
£oo - w^/('^^ - «7mw), with Eoo = 6, Wm=7.9 eV and 
7m = 51 meV. For a spherical MNP, the Green func- 
tion is computed exactlyi^. The MNP Green function 
can be understood in terms of having contributions from 
the fundamental dipolc-like plasmon mode and a reser- 
voir of higher-order plasmon modes. The dipole mode 
propagates to the far field, but the higher-order modes 
only couple in the near field. We consider an 8-nm radius 
MNP, where a QD center is at some distance h above the 
MNP surface. We assume there is a detector, at position 
td), that is 10 ^m from the surface (i.e., in the far field). 
In Fig. [2Ia), we show examples of the logarithm of 
z-projected LDOS, pdd = Pzz = lTi\[Gzz{rd,rd]Uj)] in 
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FIG. 2: (Color online) LDOS pdd at the point where the rdot 
located around a 8-nm silver MNP. (a) logarithm of LDOS at 
various locations: (top to bottom) 2 nm (red, upper), 3 nm 
(green), 4 nm (cyan), 5 nm (blue), 6 nm (magenta), 8 nm 
(black) , and 16-nm (blue dashed, lower) away from the surface 
of the MNP (from top to bottom), respectively, (b) LDOS 
pdd (blue) at Tdot situated 5 nm away from the surface of 
the MNP, and the pdo ~ \Gzz{rD,rd',iJj)\ (red), where td 
(detector position) is 10 pm above the dot position. The 
fundamental plasmon resonance of the dipole mode is ujp ~ 
2.79 eV. 



units of po, where po = fco/(67reo) is the imaginary part 
of Gji(r, r; a;) for free space. For the smallest separations 
of 2 nm (h/R = 0.25), then a QD exciton can be strongly 
coupled to the MNP resulting in vacuum Rabi splitting^; 
this strong coupling regime becomes accessible by cou- 
pling to the higher-order plasmon peaks rather than at 
the weaker dipole mode peak of the LDOS; one can see 
from Fig. [UJa) that the corresponding LDOS (red solid, 
upper) at high order plasmon modes is several times 
larger than at the dipole mode peak. For larger spa- 
tial separations between the QD and MNP surface, the 
LDOS decreased rapidly with the higher-order plasmon 
mode decreasing more rapidly than the dipole mode; the 
QD-MNP strong coupling regime is barely resolvable at 
/t, = 3 nm, and is completely lost for separation distances 
of more than h = 4 nm or h/R > 0.5. For the majority 
of our calculations we will consider h — 5 nm, though 
later we will also study the light-matter interactions at 
/i = 8 nm; for the spatial separation of /i = 5 nm, we 
have verified that a second-order Born approximation is 
valid and we will use this coupling regime to introduce a 
general master equation below. From Fig. [2][a) , one can 
see that even for h = 16 nm (h/R ~ 2), there is still an 
influence from the higher-order plasmon modes and the 
use of a single Lorentzian model would fail in general. 
In Fig. [2ljb), we show the LDOS at the h ^ 5 nm in a 
linear scale and also the magnitude of the non-local pho- 
ton "propagator," \pdD\ = \Gzz{rD,rd;ui)\. Clearly the 
reservoir function, pdd, cannot be described by a single 
Lorentzian lineshape; in contrast, the propagator, pdD, 
is mainly influenced by the dipole modeii^ii^ and is thus 



much closer to a single Lorcnzian response. 



B. Master Equation 

For the QD interactions, we consider a two-level system 
(artificial atom) in the dipole approximation, interacting 
with a general lossy and inhomogeneous structure (the 
MNP). The total Hamiltonian of the coupled system can 
be written as^i^^ 

H = h dr du}uji^{r, a;)f (r, oj) + hw^a^a^ 



/>oo 

CT+ / da;d-E(rrf,w) +H.C. 
Jo 



Hdri 



(2) 



where u^ Ig~ are the Pauli operators of the exciton 
(electron- hole pair), ujx is the resonance of the exciton, d 
is the dipole of the excitoni^, f /f '^ arc the boson field op- 
erators, and the rotating- wave approximation has been 
applied (i.e., the counter- rotating- wave term has been 
dropped). The electric-field operator (not including the 
pump field) is defined througbiii^i^i 



E(r,w) = l fdr'G{r,r';u)- J ^ei{r',u)i{r\u), 

So J V TT 



(3) 



and for convenience we have separated the 
pump Hamiltonian, defined through i^drivc 
^ (cr+e-'"^ -hcr~e*'^^), with the effective Rabi field 
fi = (Epu„jp(r£;)) • d/h. The pump field contains the 
direct pumping term plus the (dominant) scattered field 
from the MNP, 

Epump(rd,WL) = F,o{r d,ujL) + 

I dr'G(rd,r';wL)[eMNp(wL)-l]Eo(r',tJL), (4) 

where 'EiQ{vd,uji,) is the incident field operator and for 
a large driving field we can treat the Rabi field classi- 
cally (i.e., as a "c-number"). Note that the spatial in- 
tegration is carried out over the volume of the MNP; in 
this way, we recognize a plasmonic enhancement factor of 
»?? = 1 + /y^^p dv'G{Yd, r'; ujl) [eymA^L) - 1], where G 
is the total Green function of the environment including 
the presence of the MNP. Exciting near the dipole reso- 
nance with a z-polarized incident field, we estimate that 
il = Typf^o is at least one order of magnitude larger than 
fio = (Eo(rci)) • d/h. Thus the incident Rabi field is sub- 
stantially enhanced by the MNP plasmonic response22i2^. 

In a frame rotating at the laser frequency lol, we 
rewrite the above Hamiltonian as H = Hs + Hj^ + 
Hi, where the system, reservoir (or bath), and system- 



reservoir interaction terms are defined through 

Hs = h{ujx - lol)(T^cf^ + firixi(J^ + cr"), 
Hr ~ h dr dujuj{\r,uj){{r,ui), 

/•oo 

Hj = - cr+e*"^^* / <ia;d-E(rd,w) + H.c. 

"'0 



(5a) 
(5b) 

(5c) 



and then transform the Hamiltonian using H — > 
W{t)HU{t) with U{t) = exp[-i{Hs + HR)t/h], where 
the tilde denotes the interaction picture. We next ma- 
nipulate the system-reservoir interactions to derive a ME 
with the reservoir interaction included within a Born ap- 
proximation. A common choice for the ME is the time- 
convolutionless forni^i. To second order in the interac- 
tion, one has 



dm 

dt 



1 

ft2 



j dTTrR{[Hi(t), [Hi(t-r),p(t)pR 



(6) 



where p is the reduced density operator in the interac- 
tion picture, and pR = pr{Q) is the density operator of 
the photonic reservoir which is assumed to be initially 
in thermal equilibrium. Using the bath approximation, 

Trfl[f,(r,w)f/(r',w')Pfl] = H^) + l]'5(r - v')5{u: - uo') 

andTrR[f/(r,w)f,(r',w')Pfl] ="(w)<5(r-r')(5(w-cj'), we 
consider the zero-temperature bath limit (i.e., niuj) ~ 
0); which is appropriate for optical frequencies. Ex- 
ploiting the relation / dse/(s,a;)G(r,s, w)G*(s,r',a;) = 
Im[G(r, r', oj)], transforming back to the Schrodinger pic- 
ture, and carrying out the trace over the photon reser- 



= voir—, we derive the generalized ME: 



1^ = l[iJ5,p]+^*dr{jph(r)[-a+a-(-r)p + 

-hCT"(-r)pCT+]+H.C.}+/:purc(p), (7) 

where Jph(r) = /p°° dwJph(w)e'("^"'^)^, with the 
photon-reservoir spectral function given by J-p\i{uj) = 

^^^^'^ , and£purc(jO) = ■^[2criipaii-aiiaiip- 

po'iio'ii), where an = a^a" and 7' is the exciton pure 
dephasing rate. The time-dependent operators, which 
are defined through fT±(-T) = e-'-ffs^/icr^e^^^"'/'^, high- 
light that the scattering rates are pump-field dependent 
in general, since different dressed states can sample dif- 
ferent parts of the photonic LDOS^. 

For on-resonance driving (i.e., lj^ = ^x), then 



a±(r) 



= {0) 



[1 -I- cos Ht] 



'^T(O) 



cos fir] ± 



a^{0)i sin fir, which shows explicitly the formation of 
new bath-mediated scattering processes such as incoher- 
ent excitation and pure dephasing, in addition to modi- 
fied radiative decay. For numerical calculations, Eq. ([7]) 
captures non-Markovian dynamics, but for our MNP 
spectral functions, we find no evidence for non-Markovian 
behavior on the ensuing spectrum, so we can safely ex- 
tend the upper time integration on Eq. ^ to infinity. We 



subsequently obtain a useful analytic form for the ME, 



derive the following detectable spectrum, 



^ = -^ [Hs, P] + ^ {2a- pa^ - .+.-, - p.^.-) 

+ /:purc(p) + iV(r!) [a+<7- , p] + AI{n) [a+,<7'p] 

+ M*{fl) [pa'', a-] + K{n)a+pa+ + K*{n)a- pa- , 

(8) 

where the various parameters are defined as follows (see 
the Appendix for further details): 



Tin) 



N{n) = i 

M{n) = i 
K{n) = i 



d • Im[G{ojL --n) + 2G{ujl) + G{ojl + n)] ■ d 

(9a) 

d • Re[G(wL -n) + 2G{ujl) + G{ujl + ft)] ■ d 
4feo 

d • [G{ujL -n)- G{ujL + r^)] • d 

4feo 
d • [G{ujL -n)- 2G{ujl) + G{ujL +n)]-d 



(9b) 
(9c) 



4fen 



(9d) 



where G(w) = G{rd, r^; w) and we use the scattered part 
of the Green function^^. Note that in the above equations 
we have included the principal value parti^ exactly. It 
is also interesting to note that similar terms appear for 
QDs that are coupled to an acoustic phonon batb^l. 



C. Incoherent Spectrum 



SpivD^uj) ^ -\d-GirD,rd,io)\' Soiuj), (12) 

which highlights the essential role of the propagator. 
Equations ([TT |) - p2)) . together with ([7]), constitute our 
main results with which to investigate resonance fluo- 
rescence of the exciton in the vicinity of a MNP. Impor- 
tantly, the MNP reservoir function is included exactly 
and the field operators are properly quantized. An ex- 
plicit expression for the analytical spectrum is given in 
the Appendix, from which we obtain the full-width at 
half-maximum (FWHM) of the Mollow triplet center and 
sideband resonance, respectively: 



rcontcr(^) ~ y + ■^ [Jphi^L - f^) + Jph (^L + f^)] 



(13) 



rsido(^) 



J [Jphi^L -n) + 4Jph(wi,) + Jph{ujL + ^)] ■ 

(14) 

We recognize that the center line width is only affected by 
the projected LDOS (since Jph(i^) oc d-Im [G{rd,rd;uj)] ■ 
d) at the Mollow sidebands, while the sideband width 
depends on a linear combination at all three dressed res- 
onance; while such effects have been predicted before for 
atomic and dielectric syste m^^i'^° , the effects are usually 
small and to the best of our knowledge have never been 
measured nor predicted for a lossy metal environment. 



To connect to experiments on resonance fluores- 
cence, the detected spectrum will depend upon the 
position td as highlighted above for the MNP [via 
G(r£),rd)]. The spectrum is defined from S{rD,uj) = 
^(•gscatt^-j.^^ i^-)-)tE>5<=a-"(r£i, w)), where the scattering field 
due to the presence of QD is^ 



E" 



'irD,Lo) 



£0 



-G(rD,rd,a;) •d[cr (w) + cr+(a;)]. 



(10) 



For continuous wave excitation, it is common to define 
the incoherent spectrum as follows: 



^0(0;) = lim Re 



dT{{a+{t + T)a-{t))~ 



{a+{t)){a-{t)))e 



i(w£,-w)T 



where the latter term in the above equation subtracts the 
elastic (coherent) scattering from the pump field. Unfor- 
tunately, though commonly used in quantum optics, this 
expression is not valid — especially for metals — it neglects 
spectral filtering effects via light propagation. Including 
propagation and quenching in a self-consistent way, we 



III. RESULTS 
A. Asymmetric Mollow triplets 

To compute the resonance fiuorescence spectra, we as- 
sume a QD dipole moment oi d = \d\ =30 Dcbye and 
consider a pump field that excites the exciton with an 
effective pump rate, fl. Importantly, we have the dot 
in a spatial position where it necessarily feels the infiu- 
ence of the high-order plasmon modes, and consequently, 
we will show that such a system is an excellent environ- 
ment in which to study generalized reservoir coupling. 
We then solve the steady state density matrix and use 
the quantum regression theorem (or formula3^) to obtain 
the two-time correlation function and thus the spectrum. 

In Fig. [21 we display the pump-dependent Mollow 

Jx = Wp, using 



(11) spectra for the example case of uj^ 



7 



= 0.1 meV 



31 



To better clarify the role of inter- 
nal coupling effects, we show two sets of results, with 
and without the time dependence of the a^{—T) terms 
in Eq. [i.e., we set G(wl ±n) = G(wl) in Eq. dH)]. 
In Fig. [21Ja), we plot the bare spectrum with (blue solid 
line) and without (red dashed line) time dependence, and 
recognize that there is no asymmetry when the internal 
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FIG. 3: (Color online) Calculated spectra for increasing Rabi 
fields (bottom to top: Q = 20, 40, 60 meV), with ujl = co^ ~ 
ujp. (a) Spectra without propagation effects. So, with (blue 
solid) and without internal coupling effects (red dashed), (b) 
Corresponding spectra with propagation effects, Sp, which is 
the spectra that would be observed in the far field (via the 
plasmon dipole mode); to better see the sidebands, we have 
zoomed into the lower part of the spectra. 



coupling term is turned off; although there can be a small 
lamb shift, we find that this is a negligible effect. These 
findings are consistent with previous resultsii. However, 
with internal coupling, we observe a clear asymmetry in 
the Mollow triplets for sufficiently large driving fields; 
this can be explained by the complex energy (or fre- 
quency) dependence of the LDOS. For each sideband of 
the triplet, there is a one-to-one correspondence between 
the height of the sideband and the corresponding total 
transition probability between the dressed states, which 
is proportional to the LDOS at that energy [see Eqs. p^ - 
p^ ]. at the location of the QD. The large asymmetry in 
the LDOS shown in Fig.[2]thercfore manifests itself in the 
asymmetric strength of the Mollow triplets for sufficiently 
large drive fields. For the case without internal coupling, 
assuming large drives, then only the LDOS at the fre- 
quency of the pump field matters and the corresponding 
Mollow triplet is symmetric. We also observe significant 
narrowing of the resonances, consistent with Eqs. ((13])- 
P^ ]. In Fig.|31Jb), we observe similar features for Sp{ijj), 
but now the spectra arc further reshaped due to photon 
propagation and quenching effects; specifically, without 
propagation effects the MTs are symmetric without inter- 
nal coupling, and asymmetric with IC; with propagation 
effects, both are asymmetric, but in opposite senses. 

Next, we investigate the case when the pump field 
is resonant with the exciton, but off resonant with the 
dipole mode of the MNP, and show that the spectral 
filtering effect can be used to selectively enhance the fea- 
tures of the Mollow sidebands. In Fig. SI we show the 
high-field solution when the pump field is now 60 meV 
blucshiftcd with respect to the plasmon dipole mode 
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FIG. 4: (Color onhne) As in Fig. H using fi = 60 meV, 
but now the pump frequency is 60 meV above the fundamen- 
tal plasmon mode: ul = oOp + 60 meV. (a) Spectra without 
propagation effects, 5*0, with (blue solid) and without internal 
coupling effects (red dashed), (b) Corresponding spectra with 
propagation effects. Here the fundamental plasmon mode is 
seen to relatively enhance the lower Mollow sideband. 
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60 meV). Here the photon propagator 



has changed the spectrum of the scattering field dramati- 
cally. Due to the fact that the left sideband predominates 
over both other peaks, it should be more accessible in an 
experiment. Figure IH^a) shows the following features: 
(i) the induced- asymmetry from internal coupling is now 
minor since the LDOS at ujl i 60 meV is similar, and 
(ii) there is significant spectral broadening when internal 
coupling is included; this broadening effect is anticipated 
from our theory as the Mollow sideband resonances now 
sample a larger LDOS [see Fig.[T]and Eqs. ((131) -(HI)]. We 
stress that such effects are not possible in a Lorentzian- 
decay medium (e.g., a single mode cavity). 



B. Position dependence of the Mollow triplet 

Due to the fact that high-order plasmons are strongly 
confined near the surface of the MNP, one may expect 
that as the distance between the QD and MNP is in- 
creased, then a simple single Lorentzian may be valid. In 
addition, it is useful to know by how much the Mollow 
triplet features change if one moves the spatial position 
of the QDs by a few nm. As is shown in Fig. [2ja), even if 
the distance between the QD and the surface of the MNP 
is as large as the radius of MNP, the higher-order plas- 
mon modes still have a significant impact and the failure 
of a single Lorenzian model is to be expected. In fact, 
we find the photonic resevoir of QD induced by the MNP 
will always display some non-Loretzian characteristics. 

Figure [5] presents the Mollow spectra with the pump 
frequency lol = ^x = Wp. Figure ^a) shows the bare 
spectrum with (blue, solid line) and without (red, dashed 
line) internal coupling, and clearly there is still an asym- 
metry between the Mollow sidebands; however, since the 
LDOS values are reduced the magnitude of the sidebands 
arc suppressed and the radiative decay rates are much 
smaller (less broadening). In Fig. [5l^b), we show the de- 
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FIG. 5: (Color online) Calculated spectra for the QD located 
8-nin away from the surface of the MNP with increasing Rabi 
fields (bottom to top: Q, = 20,40,60 meV), and ujl = uJx = 
ujp. (a) Spectra without propagation effects. So, with (blue 
solid) and without internal coupling effects (red dashed), (b) 
Corresponding spectra with propagation effects, Sp. 



tectable spectrum in the far field, which again shows sig- 
nificant reshaping of the spectrum due to photon prop- 
agation and quenching effects. However, as the spatial 
distance between the QD and MNP is further increased, 
the contribution of higher-order plasmon modes becomes 
smaller and smaller; eventually we find that when the QD 
is placed 16 nm (2i?) away from the MNP surface [blue 
dashed in Fig. [H^a)], then the magnitude of the LDOS 
can be reasonably described by the dipole mode only as 
the LDOS contribution from the higher-order plasmon 
modes is about one order of magnitude smaller than the 
dipole mode peak LDOS; in this larger separation regime, 
the MNP can be effectively described by the dipole model 



as also discussed elsewhere^*^. 



IV. CONCLUSIONS 

We have introduced a general ME approach for 
modelling quantum light-matter interactions for a driven 
atom or QD in the vicinity of a metallic nanoparticle. 
The exact reservoir function and propagator are obtained 
from Green function theory and used directly in the ME 
formalism. We used this approach to model the MoUow 
spectrum as a function of drive strength, and demon- 
strated that clear spectral asymmetries and non-trivial 
linewidth variations can be seen for suitably large drives. 
We also investigated several different pump excitation 
frequencies and QD positions and find rich coupling 
behaviour. While our master equation formalism is 
useful for exploring regimes of quantum nanoplasmonics, 
the techniques are general and can, assuming the 
validity of the second-order Born approximation, model 
the spectrum from a driven QD in any general photonic 
environment, including hybrid metal-photonic-crystal 
systems^ and metamaterial o^^i'^^ . 
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Appendix: Derivation of the Incoherent Spectrum and Full- Width at Half Maximum of the Mollow Triplet 

Resonances 

Using the master equation, Eq. ([8|), we derive the following Bloch equations 



d(a+) 
dt 

d{<J-) 

dt 
dia^ 

dt 



l^mttJ^ + N{n}^ {a+) + K*{n){a-) - ^^(a^) - A/*(0), 
(2A/(f7) + in){a+) - (2Ar(0) - in){a-) - r(f7)((CT^} + 1), 



(A.l) 
(A.2) 
(A.3) 



with parameters T{n) defined through Eq. (Pa)) - (j9dp in the main text of the manuscript. These Bloch equations can 
be solved exactly, which is the approach we have used in the main text above, with no approximation. However, it is 
useful to look at certain limit of the corresponding solution for the spectrum, 6*0(0;). This allows us to make a clearer 
connection to the underlying physics of the resulting spectral linewidths. Neglecting the terms from the real part 
of the Green function (since these only cause spectral shifts through the Lamb- and Stark- shifts), the steady-state 



solutions of the Bloch equations are as follows: 



(a+) 



+ \ -,■ 



^ /r(oi+y _ ^^(^) _ ^M^ 






r(o) 



Mrin) 



2 ( iiiiii^ - X, 



r(o)+7' 



,(fi)) ( 



r(n)+7' 

2 



.fi 



Ar(ii) + pfoN 



r(n) 



2 



Xr(f7) ' 



Mrin)+t-] (a+)ss 



1, 



(A.4) 

(A.5) 
(A.6) 



where Kr{^) ~ Rc[K{il)] and A/r(ri) = Re[M(r2)]. The incoherent spectrum, without propagation effects (see main 
text), is given by 



So{uj) = lim Re 



dr((a+(f + r)<7-(i))-(a+(i))(a-(i)))e^(---")^ 



Re^ 



' d{0) X -^KAn)Su.-2MAn)n-m^+2KAn)rm ^ ^^^^^ _ 3^^^ ^ jr(17)]/(0) ' 



2Suj(Kr{n) 



3r(0)+7' 



+ 2i 



-172 ^ 5^2 



X^(0) 



r(o)+7 



-) r(fi) 



(A.7) 



with d(0) = /(O) +g(0), /(O) = (dV+^cr"), 5(0) = {5a+6a+), and /i(0) = {Sa+Sa''), respectively. These terms are 
obtained from the steady-state Bloch-cquation solutions, via 

/(O) =^(1 + ('T^)s. - 2(a-),,(<7+)33), 5(0) = -(^+)L /i(O) = -(a+),3(l + {a%^). 

In the strong pump limit (i.e., fl > 7', \T{n)\, |iV(0)|, |M(f})|, |i^(rj)|), the full-width at half maximum (FWHM) 
values of the MoUow triplet resonance are obtained from the imaginary parts of 



r(f}) + 7' 



iSu} — A'r(r2) 



2i 



Slu^ - n^ 



Krin)- 



T{n) + j' 



nn) 



-2Slu 



Krin) 



The corresponding roots are easily obtained: 



Sluo 



6uj+ = ± 






Krin) , 



4 [02 + r{n)B{n)] - [r(o) + B{n)Y 



i[T{n) + B{n)], 



3r(o) + 7' 



(A.8) 

(A.9) 
(A.IO) 



with B{Q,) — ^ ^^^ + AV(0). In the strong field limit, the real parts of the roots correspond to ujl,^l ± ^, at these 
MoUow triplet resonance, and the FWHM of the spectral linewidths are 



^ center t,^ ' j 
rsido(f^) = TI" 



7' , Jph(t^L - f^) + -'"ph(wL + 0) 

h TT , 

2 2 

jph(^L - 0) + 4jph(^L) + jph(^L + ^) y 

4 4' 



(A.ll) 
(A.12) 



which explicitly show the role of the three LDOS values at the dressed-state resonance. 

It is also useful to compare the above bare-state approach with an approximate stressed-state approach in the 
secular approximation, which has been used before in the context of coupling to generalized reservoirs^Si^. The 
Mollow triplet can then be explained from the energy level scheme in the dressed-state picture, where the MoUow 
central peak is due to the evolution of {(J^^^ss (*) ) ' ^^'^ ^^^ sidebands are related to the relaxation of the dipole operators 
((T^^ggg(i)). The dressed-state operators are related to the bare state operator through the following transformations. 



Mre 



K^z 



{a^ +a^ —fj^) and a^ 



a~^ +a . Using bare state operators, and adopting the secular approximatio: 



,29 



we obtain the following Bloch equations, 






dt 



dt 



Jph(wL - ri) - 2jph(u;L) + jph(wL + ri) jph(wL - 17) - ^ph(wL + ri) 7' _j_, 

7^ ^^ (o-+)+7r — (cr ), 

/ Jph(wL - f7) - Jph(wL + ri) - , . Jph(wL - 11) - Jph(wL +17) . _ 

= (7r til) {a } + [TT hii2j((T )- 

-'"ph(wL - 17) + 2Jph(wL) + JphCl^L + 17) ., z. , ,s 
TT 7^ ((ct ) + 1). 



(A.13) 



(A.14) 



Thus the time evolution of the average values for the dressed state operators are given by 



dt 

dt 



I z .TriJphCll^L - 17) + Jph(t^L + 17)] +7' , Jph(wL - 17) - Jph(wL +17) 
(C^drcss) ^ TT- 



(crr ) I ±il7 

^ dress/ ' 



2 2 

TT [Jph(wL - 17) + 4 Jphlwi) + Jph(wL + 17)] + 7' 



(^drcss) 



TT [JphCl^L - 17) + Jph(wL + 17)] + 7' 



_ Jph(wL - 17) + 2Jph(wL) + Jph(wL + 17) 



(O-drcss/^ 



^ph(wL - 17) - Jph(l^L + 17) 



(A.15) 



(A.16) 



from which we obtain the spectral linewidths 



,f^. _ 7' , >/ph(wL - 17) + Jph(WL + 17) 
ccntori"j — -^ + TI" ^ , 

7' , ^/phCl^L - 17) + 4Jph(wL) + Jph(wL + 17) 



rcoh(17) = rsidc(17) = — + TT 



(A.17) 
(A.18) 



As expected, these are in agreement with the more exact bare state approach when analyzed in the MoUow limit 
(which is similar to making the secular approximation). However, the advantage of Eq. ^ is that it can be applied 
for all values of the pump field, so the secular approximatio n is not needed. 
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